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Abstract 

The inner structure of the attractor appearing when the Varley-Gradwell- 
Hassell population model bifurcates from regular to chaotic behaviour is 
studied. By algebraic and geometric arguments the coexistence of a con- 
tinuum of neutrally stable limit cycles with different periods in the attractor 
is explained. 
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1. Introduction 

Since the first discoveries of chaotic behavior in discrete maps, there has 
been great interest in the description of how the dynamics of a map changes 
from regular to chaotic, and viceversa, when a parameter of the map is 
varied. The so-called routes to chaos have been widely studied for smooth 
maps where well-known phenomena such as period-doubling, saddle-node 
bifurcations or intermittency appear jH 121 Ej • The study of piecewise defined 
maps, which are useful in the mathematical description of mechanical systems 
with friction or collisions and of electrical circuits with switching components, 
has introduced new phenomenology [H [5] . 

Here we consider a one- dimensional piecewise smooth discontinuous map 
originally proposed in studies on population ecology [SIE], the Varley-Gradwell- 
Hassell map (henceforth referred to as VGH map). It is a three-parameter 
map composed of a linear part and a power-law decreasing piece. It may 
present both regular and chaotic behavior, and for a critical value of one of 
its parameters the system undergoes an abrupt order-to-chaos transition. In 
[8] it was shown that at this transition point, the system has a continuum of 
neutrally stable limit cycles, i.e. cycles with multiplier unity. That situation 
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appears also in two-dimensional maps and in a recent work [H] we have named 
such bifurcations bisecting bifurcations for a reason that will become clear 
in the following sections. The interesting characteristic of the bisecting bi- 
furcation in the VGH map is that neutrally stable limit cycles with different 
periods coexist at the transition point. Their periods and their distribu- 
tion in phase space depend on the location of the discontinuity of the map. 
Our previous results are here sharpened by invoking algebraic and geometric 
arguments. 

The structure of the paper is as follows. In Section [2] we introduce the 
VGH map and its bisecting bifurcation. In Section |3] we explain a way to 
study algebraically the existence and location of this type of bifurcations. In 
Section [4], using the cobweb diagram of the map we provide a geometrical ex- 
planation of the bisecting bifurcation as well as the reason for the coexistence 
of different periods. In Section |5] we detail the structure of the attractor and 
of the basins of attraction at the bifurcation point. And finally, in Section [6] 
we discuss the main results presented in the paper. 

2. The VGH map 

The VGH map is defined in [0, oo) and reads 

^^^^ ~ I rx^~^, X > c, ^^"^ 

where r,b > 1 and c > 0. When c = 1 the system is continuous, being 
discontinuous otherwise. The VGH map presents both regular and chaotic 
behavior and its dynamics exhibits a number of phenomena including interior 
crises and crisis-induced intermittency [8]. 

A customary construction of the bifurcation diagram by varying b, for 
fixed r and c, shows an order-to-chaos bifurcation at 6 = 2. For b < 2 the 
system is regular with a single periodic attractor whose period depends on 
c. When b > 2 the system is chaotic without further regular window. This is 
clear from the bifurcation diagrams in Figure [T] where three different values 
of the parameter c are considered. The top and bottom panels of the figure 
stand for cases where the map is discontinuous (c = 0.8, 1.2) and the middle 
panel shows the bifurcation diagram in the continuous case (c = 1). At the 
bifurcation point 6 = 2 all these diagrams exhibit a vertical segment which 
constitutes the focus of this work. 



2 



Detailed numerical experiments, buttressed by the arguments in the next 
sections, demonstrate that any initial condition end up in a neutrally stable 
limit cycle. There is precisely an infinity of them which form the vertical 
segment in the diagrams of Figure [1} It is worthwhile to point out that the 
numerical detection of this continuum of attractors which mediates the bifur- 
cation presents three main difficulties. First of all it exists only for the precise 
numerical value 6 = 2. Then, an accurate enough machine-representation of 
the parameter values must be used in the sampling. Secondly, as the num- 
ber of neutrally stable periodic orbits that exist for 6 = 2 is infinite, a large 
enough number of initial conditions must be used to generate the vertical 
segment in the diagram. Finally, the whole attractor can easily be mistaken 
for a single chaotic attractor unless properly scrutinized. The algebraic and 
geometric analyses in the next sections are of great help in the correct iden- 
tification of the vertical segment in the bifurcation diagram. 



3. Localization of the bifurcation in parameter space 

In this section we prove that for the VGH map the bisecting bifurcation 
takes place only at 6 = 2. We start by introducing a change of variables that 
permits us to exactly linearize the VGH map 

z = 21og(a;)/log(r), ^ = 2 log(c)/ log(r). (2) 

In terms of z and the VGH map can be cast in the following linear form 

This linearized version allows the following algebraic approach. In the new 
form L{z\ the elements {zi, ^2, • • • , ^t} of a cycle of period T satisfy 

Zn-\ + 2, Zn-\ < ^, 

(1 - \))Zn-\ + 2 Zn-\ > n = 2, . . . , T, 

_ r 2;t + 2, ZT<i, .^x 

1 (l-6)zT + 2, ZT>i. ^ ^ 

These expressions can be cast in the form 

Zl + axZT = 2, 
Zn + an-iZn-^i = 2, n = 2, 3, . . . , T, (5) 
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with coefficients 

(6-1). 

The hnear system of equations ([s]) can be expressed in matrix form 

T 

Y,AjZ, = 2, z = l,2,...,T, (7) 

i=i 

where the matrix A is 

/ ai 1 \ 

aa 1 ... 
. 1 ••• 



ar-i 1 

\1 ar J 

In order to find an attractor consisting of infinitely many hmit cycles of 
period T we need the system of equations ([s]) to be compatible and indeter- 
minate. A necessary condition is therefore det(y4) = 0. The determinant of 
the matrix A is 

T 

det(A) = + {-lf+\ (8) 

fc=i 

If the cycle has d of its elements satisfying Zi > ^ this gives 

det{A) = {-lf-\b - lY + (9) 

which vanishes for 6 = and, if d is even, also for 6 = 2. Since we are only 
interested in 6 > 1, the only possible solution is then 6 = 2. This is precisely 
the value found from the bifurcation diagrams in Section |2} It is important 
to notice that this value of 6 and equation ^ imply neutral stability for any 
possible cycle. Of course, the argument in this section does not complete the 
analysis of the vertical segment in those diagrams. It only excludes any other 
value of 6 as candidate. In the next Section we will explicitly show that the 
infinite set of periodic orbits actually exist. 
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4. Geometrical explanation of the bisecting bifurcation 



Next we develop an heuristic explanation for the emergence of a bisect- 
ing bifurcation based on the cobweb diagram of a general one-dimensional 
piecewise defined map f{x; p) with parameter p. Notice that an infinite set of 
regular trajectories of period n will be observed for a critical value pc of the 
parameter when the nth iterate of the map /["^ has a piece which is co-linear 
with the bisectrix. Hence, the name we are using for these bifurcations. More 
explicitly, at the bifurcation, /'"l will be of the form 

= I xe {Xi,Xr), (10) 

An instance of this phenomenon in the VGH map is given by its second 
iterate which reads 



r'^'^x^'^, ^ < X < c, 
r^x^~^, x>c,x>{-)^, 



v^'\x) = { :^^Zir_,^ (11) 



When 6 = 2 the third piece becomes x and we find a bisecting bifurcation 
mediated by an infinite set of neutrally stable period-2 limit cycles provided 
c < a/t. This is illustrated in Figure [2] where three cobweb plots of V^^ix) 
for values h = 1.6, 2.0, 2.4 and c = 1 are shown. In the three cases the same 
three initial conditions have been evaluated. When b = 1.6, all trajectories 
converge to the stable fixed point. In the critical case, when 6 = 2, each 
initial condition evolves to a different neutrally stable fixed point of V^'^^{x). 
For b = 2.4 the trajectories are chaotic. 

In principle, several iterates of / can have pieces co-linear with the bisec- 
trix simultaneously for the same value of pc- In such cases, infinite sets of 
cycles of different periods will coexist when p = p^. This is the case for the 
VGH map. When 6 = 2 and c > 1 the neutrally stable limit cycles can show 
different periods. This is illustrated in Figure [3] where the evolution in the 
cobweb diagram and the trajectories of two different initial conditions are 
shown for 6 = 2, r = 4 and c = a/2. One of the trajectories enters a period-2 
limit cycle while the other evolves to a period-4 limit cycle. The distribution 
of the limit cycles of different periods in phase space is detailed in the next 
Section. 
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5. Attractors and basins of attraction for b = 2 

In this Section we study at deptli tlie structure of tlie attracting segment 
of the VGH system at the critical point 6 = 2. To facilitate the analysis we 
express the VHG map for 6 = 2 in terms of the new variable w = z — 1 as 

W{w) = l ^ + 2, ^<^-|' (12) 
^ ^ [ -w, w > ^ -1. ^ ^ 

One advantage of this form is that the segment of limit cycles, A = [— 1^| — 
1, 1^1 + 1], is symmetric in phase space with respect to w = 0. 

5.1. Structure of the attractor 

We start our study by distinguishing the cases of positive and negative ^. 

5.1.1. ^ < 

This case presents the simplest dynamics. The segment A is composed 
of infinite period-2 limit cycles around w = 0, which is a fixed point. 

5.1.2. ^ > 

The dynamics is more complicated for positive ^. All integer initial con- 
ditions Wo G Z lead to limit cycles with integer elements. In particular, if 
the initial condition is even (resp. odd), the final limit cycle will have as its 
elements all even (resp. odd) integers inside the segment A. The periods of 
these cycles, which depend on ^, are detailed in Table [T} If ^ Z we need 
to study the cases ^ € N and ^ ^ N separately: 

• e ^ N 

For positive non-integer ^ we have two infinite sets of limit cycles of 
periods T = 2(n + 1) and T = 2(n -|- 2) with n = [^J (where [-J stands 
for the fioor function). These cycles spread over ^ in a rather peculiar 
way. To make it clear we find useful to consider the finite set of points 

B = Af]{w = {-l)''^+{-iy{2k + l), k = 0...n, a, (3 = 0,1} 

(13) 

whose elements, when written in increasing order, we denote hj w\ i = 
1, 2, . . .. Then, B punctuates a partition of the interval A in subinter- 
vals. Points in the same subinterval belong to cycles with the same 
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period. Points in contiguous subintervals belong to cycles with differ- 
ent period. 

For instance, if < ^ < 1 then n = and the frontiers of the subinter- 
vals are given by 

B = {w\ w\ w\ w'} = {-e - 1, e - 1, -e + 1, e + 1}. (14) 

Thus, in this case three subintervals exist inside the attractor. In Table 
[T] we have detailed the period of the limit cycle for initial conditions 
inside the attracting segment A. 

When ^ is a natural number every limit cycles has period T = 2(,^ + 1). 

All this information is contained in Figure |4] where we have plotted the 
structure of the attracting segment at 6 = 2 as a function of ^ with ^ G [—1, 3]. 
The different tones of gray stand for different periods. As an example, a 
particular ^* < 1 has been chosen to illustrate the position of the frontiers 



given in (14). The discontinuous horizontal lines stand for cycles with initial 



conditions Wq G Z. 

5.2. Basins of attraction 

Given the limit cycle to which a point wq tends, its basin of attraction 
can be written as 

= {w e M.\w = [wo\ +2k±d,k eZ}, d = wq - [wo\ . (15) 

This structure of the basins of attraction is reflected in Figure |5] for two dif- 
ferent values of ^. In this figure for each initial condition wq G [—3, 3] we 
plot the cycle in A which traps it. The figure suggests a periodic structure 
in the horizontal direction. It reflects the partition of phase space into equiv- 



alence classes established by (15). More specifically, all even, odd and odd 
half-integer initial conditions constitute three equivalence classes by them- 
selves. Any other initial condition generates its equivalence class by repeat- 
edly adding alternatively 2d and 2(1 — li). 



The combination of Table 1 and (15) allows to determine the period of 



the cycle to which an arbitrary point wq tends. For the sake of illustration. 
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consider the case ^ = 1.5 and the initial condition wq = 2.7. In this case 
n = [^\ = 1 and the frontiers inside the segment A are given by 



B = w^ U7^ w\ w^} = {-2.5, -1.5, -0.5, 0.5, 1.5, 2.5}. (16) 

Since the chosen initial condition falls outside the attracting segment we will 
make use of (15) to determine another initial condition w'q G B^^ leading to 
the same final limit cycle. Choosing A; = we can readily find 

w'q = [wo\ +2k-d=1.3e {w\ w^). (17) 

If we now take into account the classification detailed in Table [T] we can 
conclude that both wo and w'q will enter a limit cycle of period 4. It can be 
checked that this is in fact the case in the upper panel of Figure |5| 



6. Discussion of the results 

This work focuses on the coexistence of different periods in the set of 
neutrally stable limit cycles that mediates the bisecting bifurcation of the 
VGH map. We have proved algebraically that it occurs only for 6 = 2. In 
principle, the same method might be used to determine the existence and 
position in parameter space of these bifurcations for any piecewise linear 
map. 

The geometric explanation complements the algebraic approach and pro- 
vides us with an intuitive way of interpreting the bifurcation in terms of the 
cobweb diagram. In particular, it allows us to understand the coexistence 
of different periods as the simultaneous presence of pieces co-linear with the 
bisectrix in different iterates of the map. 

The structure of the attractor in 6 = 2 has been described in detail. Given 
a certain value of the discontinuity parameter (either c or ^), the periods of 
the limit cycles and their location in phase space can readily be determined. 
The structure of the basins of attraction of the limit cycles has been resolved 
as well, allowing us to associate each initial condition to its final limit cycle. 

The described bisecting bifurcation is accompanied by border collisions 
[H El [in] as can be seen in Figure |2} In the example described in that figure, 
two unstable fixed points of V^'^^{x) (for 6 > 2) collide with discontinuities 
of the map when 6 = 2 and cease to exist for 6 < 2. An attempt to clas- 
sify the phenomenology of border-collision bifurcations in one-dimensional 
discontinuous maps can be found in [TT]. The classification is based on the 
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linearization of the map around the colhsion point both in phase space and 
parameter space. However, the bisecting bifurcation presented here corre- 
sponds to one of the critical cases explicitly excluded from this classification. 

In this study we have restricted our attention to the VGH map but what 
we have called bisecting bifurcations are present in other continuous and 
discontinuous piecewise smooth maps. However, to the best of our knowledge, 
they have very often gone unnoticed in the literature. In the continuous case, 
bisecting bifurcations can be observed in maps such as the skew tent map [12] 
or the map describing the dynamics of the boost converter [13] . In particular, 
in [13] a continuous piecewise smooth map introduced as a model of economic 
growth [15] is studied and the values of the parameters for which a bisecting 
bifurcation takes place are identified. Discontinuous maps candidates to show 
bisecting bifurcations can be found in [161 HZ] ■ Our numerical experiments 
have shown this is in fact the case. Moreover, the map studied in |l16J shows 
coexistence of different periods in the set of neutrally stable limit cycles with 
a structure very similar to the one described in this paper for the VGH map. 
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Table 1: Period of the trajectories according to the initial point wq E [~\^\ — 1, |CI + 1] 
with 6 = 2 and ^ > ^ N). The description of the case wq E Z is also valid for ^ G N. 



Initial point 


Period 


Condition 


Wq G {w\w'^~^^) and wq ^ 


2(n + l) 
2{n + 2) 


even i 
odd i 


Wq = and Wq ^ Z 


2{n + 2) 


Vi 


Wq E Z 


^=L(^ + 1)/2J 2N + 1 
M = [^/2J 2(M+ 1) 


even wq 
odd Wq 




Figure 1: Bifurcation diagrams of the VGH map V{x) with r — 2 and three different 
values of c. 



11 




Figure 2: A sequence of three cobweb plots of T^I^l (x) illustrating the bisecting bifurcation 
from order to chaos at 6 = 2. Here c = l,r = 2. The presence of the segment of fixed 
points is apparent in the case 6 = 2. 




Figure 3: Cobweb diagram (left panel) and trajectories (right panel) of V{x) for 6 = 2, 
c = v^, r = 4 and two different initial conditions (a;o = 0.3,0.6). It is clearly seen that 
each initial condition enters a neutrally stable limit cycle of different period. 



12 




Figure 4: Structure of the attracting segment {b = 2) in phase space for different values 
of ^. The diffcrcint tones of gray code cycles with different periods. The discontinuous 
horizontal lines stand for individual cycles with integer elements of periods 1, 2, 3 and 
4 (dashed, dotted, dash-dotted and dashed respectively) embedded in the continuum of 
cycles of the attracting segment. 




Figure 5: Structure of the basins of attraction for two different values of the discontinuity 
location parameter ^. The graphs show the final attractor for each wq. For ^ = 1.5 (upper 
panel) infinite limit cycles of periods 4 and 6 exist as well as the period-2 cycle {-1,1} and 
the period-3 cycle {-2,0,2} can be seen. For ^ = 2.0 (lower panel) infinite limit cycles of 
period-6 arc present filling the space between the points of the period-3 cycle {-2,0,2} and 
the period-4 cycle {-3,-1,1,3}. 
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